Re-entrant ferromagnetism in a generic class of diluted magnetic semiconductors 
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Considering a general situation where a semiconductor is doped by magnetic impurities leading 
to a carrier-induced ferromagnetic exchange coupling between the impurity moments, we show 
theoretically the possible generic existence of three ferromagnetic transition temperatures, Ti > T2 > 
T3, with two distinct ferromagnetic regimes existing for T\ > T > T2 and T < T3. Such an intriguing 
re-entrant ferromagnetism, with a paramagnetic phase {T2 > T > T3) between two ferromagnetic 
phases, arises from a subtle competition between indirect exchange induced by thermally activated 
carriers in an otherwise empty conduction band versus the exchange coupling existing in the impurity 
band due to the bound carriers themselves. We comment on the possibility of observing such a re- 
entrance phenomenon in diluted magnetic semiconductors and magnetic oxides. 

PACS numbers: 75.10.-b ,75.50.Pp , 75.30.Hx 



The large class of materials comprised of diluted mag- 
netic semiconductors and magnetic oxides (DMS) has 
been widely studied in the recent years for their potential 
in spintronic applications as well as for the fundamental 
physics of carrier-mediated ferromagnetism in semicon- 
ductors. They have ferromagnetic (FM) critical temper- 
atures Tc ranging from below (e.g. (Ga,Mn)As (ij) to 
above room temperature (e.g. doped magnetic oxides as 
Ti02 0), and some others never show long range FM 
order (e.g. most Mn-doped II- VI semiconductor alloys 
are spin glasses) Q. FM is usually ascribed to carrier- 
mediated mechanisms and depends on many different 
parameters (e.g. carrier density ric, magnetic impurity 
density Ui, magnetic coupling between the ion moment 
and the electron (or hole) spins J, and details of disor- 
der) which vary greatly from system to system and some- 
times from sample to sample. When the concentration of 
magnetic impurities is larger than a few percent, direct 
exchange (i.e. not carrier-mediated) can also be signifi- 
cant. This exchange is antiferromagnetic in many cases, 
like the Mn-doped III-V [3] and II- VI semiconductors Q , 
and ferromagnetic in others, like Co-doped Ti02 0- 

The accepted effective theoretical models 0, H, @, [13, 
ll| for FM in DMS can be roughly divided into two broad 
categories depending on whether the carriers mediating 
the FM interaction between the magnetic impurities are 
itinerant free carriers (i.e. holes in the valence band or 
electrons in the conduction band) or localized (or bound) 
carriers (in an impurity band, for example). For the itin- 
erant free carrier case, e.g. Gai_a;Mna;As in the opti- 
mally doped (x Ri 0.05) situation with the nominally 
high Tc 100 — 170 K), the effective magnetic inter- 
action producing FM is thought to be of the Rudernian- 
Kittel-Kasuya-Yosida (RKKY)-type [3, i, [i, [S, [ill [11 
leading to a mean-field FM transition by virtue of there 
being many more effective impurity magnetic moments 
than free carriers. Such RKKY carrier-mediated mean- 
field FM seems to describe well 0, [13, [m the optimally 



doped (Ga,Mn)As regime. For the localized case, where 
disorder should play an important role, the FM transi- 
tion is thought to be caused by the temperature-driven 
magnetic percolation transition of bound magnetic po- 
larons (BMP) [l^, objects where one localized carrier is 
magnetically strongly correlated with a few neighboring 
magnetic impurities through the local exchange interac- 
tion. An example of such BMP percolation induced FM 
is thought to be the localized impurity band Gei-xMn^ 
DMS system [li|. Typically, the RKKY (BMP percola- 
tion) ferromagnetic Tc is relatively high (low). 

In this Letter, we argue that it is possible, perhaps 
even likely in some impurity band systems, for these two 
mechanisms to operate together in (at least some) DMS 
materials, leading to an intriguing situation where the 
RKKY mechanism, mediated by thermally excited free 
carriers, dominates at high T, whereas the low T regime 
(where thermal carrier activation is exponentially sup- 
pressed) is dominated by the polaron percolation mech- 
anism in the impurity band. By theoretically analyz- 
ing such an interesting multi-mechanism DMS situation 
using simple physical models, we show that it may be 
generically possible for DMS materials to have re-entrant 
FM with a paramagnetic state sandwiched between the 
high-temperature 'activated' RKKY-type and the low- 
temperature BMP percolation-type FM. 

The system we consider is a localized impurity band 
system which is insulating at T = 0. In this scenario, the 
density of band carriers is a function of temperature T 
and activation energy A, 



nc{T, A) = Uco exp (- A/ZcbT) 



(1) 



As a consequence, the RKKY mechanism is inefficient at 
low T (typically for T < A) when the free carrier den- 
sity is too low, while the activated carriers can mediate 
FM at higher temperatures. Therefore, this 'activated' 
RKKY mechanism leads, by itself, to a FM phase be- 
tween two disordered phases at high T > Ti and low 
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FIG. 1: (a) Impurity magnetization {Sz)/S corresponding to 
the 'activated' RKKY model, Eq. ([2]), for different values of 
the carrier activation energy A. The parameters used are 
rico = 3 X 10^° cm-^ n, = 3 x 10^^ cm-^ J = 1 eV. (b) The 
two critical temperatures Ti and T2 of the 'activated' RKKY 



model, Eq. ([5}. T and A are in units of Tj( 



Ti(A = 0). 



T < T2 temperatures. At low T, other mechanisms not 
requiring free carriers (e.g. magnetic polaron percola- 
tion) can come into play producing a further magnetic 
transition with critical temperature T3. When T3 < T2 
there are two distinct ferromagnetic phases and the sys- 
tem exhibits re-entrant FM. 

The Hamiltonian of the exchange interaction between 
magnetic impurities and the carriers is H ^ J2i Jo,a ' 
s(Ri) where J is the local exchange coupling between 
the impurity spin Si located at Ri and the carrier spin 
density s(r), and Uq is the unit cell volume. 

As an illustration, we first discuss the hypothetical case 
of an intrinsic semiconductor with magnetic impurities 
which do not contribute carriers (therefore, there exists 
no impurity band where BMPs could form). In this case, 
the RKKY mechanism can only be mediated by carri- 
ers thermally activated from the valence to the conduc- 
tion band. Another mechanism that could play a role in 
this intrinsic case is the Bloembergen-Rowland mecha- 
nism [1^ mediated by virtual electron excitations across 
the band gap. This mechanism was considered in the con- 
text of II- VI DMS [3j] and was finally dismissed as being 
too weak compared to superexchange Therefore, we 
do not further include this mechanism in our discussion. 

We now consider this 'activated' RKKY mechanism 
within a mean-field theory in the limit of nondegener- 
ate carriers. In the activated carrier scenario we are 
studying, the Fermi energy is smaller than the tempera- 



ture, and hence the limit of nondegenerate carriers is ap- 
propriate. Within mean-field theory, the impurity spins 
act upon the carrier spins as an effective magnetic field 
oc Ja^ni{Sz) while the carrier spins act upon the impu- 
rity spins with an affective field oc J a^nc{s z) ■ As a result, 
the magnetization of the magnetic impurities {Sz)/S is 
calculated self-consistently giving 



(Sz) 

s 



Be 



Jaln,iT,A) ( Jaln,{S,) 



(2) 



where Bs{y) is the standard Brillouin function. The mag- 
netization is plotted in Fig.[l](a) for three different values 
of A. The A = curve has a non-standard concave shape 
as expected for a low density of carriers . The critical 
temperatures are given by the points at which the mag- 
netic susceptibility diverges. The susceptibility is essen- 
tially that of the magnetic impurities, as usually Uc << 
rii due to carrier compensation, and in the paramag- 
netic regions is given by x(r) = {gifJ-siSz)) /OB. For 
A = 0, there is only one critical temperature given by the 
standard result T^^^^ = \ j a%^/n;:^^^ [S + l)(s -\- 1). 
For A 7^ 0, the model gives two critical temperatures T\ 
and T21 shown in Fig. [1] (b): at low T, ric gets expo- 
nentially suppressed and the localized moments are in- 
dependent of each other; as T increases, the band gets 
populated and the carrier-mediated FM kicks in at Ti-, 
which increases with A. At even higher T, thermal dis- 
order produces the standard ferro-to-paramagnetic tran- 
sition with critical temperature T\ , which decreases as A 
increases, due to the reduction of carrier density in the 
band. The Ti and T2 curves meet at A « 0.73 T^^^^"^ 
and, for A > 0.73 T^^^^, the density of carriers is al- 
ways too low to mediate FM so the 'activated' RKKY 
model gives paramagnetism at all T. The curves for dif- 
ferent parameters scale with T^^^ = ri(A = 0). 

We consider now the more realistic case in which there 
is an impurity band. The carriers in the impurity band 
can come from magnetic impurities acting as dopants, 
like in the III-V semiconductors, or from other dopants, 
like the oxygen vacancies which act as shallow donors in 
the magnetic oxides. We assume that, at T = the con- 
duction (or valence) band is empty so all the free carriers 
in the band are produced by thermal activation from the 
impurity band. At low T, the carriers in the impurity 
band can be strongly localized and mediate FM through 
the formation of BMPs which grow as T lowers, finally 
overlapping in clusters and producing a FM transition 
at percolation In our model of thermally activated 
carriers, the density of carriers in the impurity band is 
n*(T, A) = Uca ~ = tIco [1 - exp(-A/fcBr)], where 
Tif, is the density of carriers in the valence or conduction 
band. A is the smallest energy gap for electron excita- 
tions, which is given by the activation energy of a carrier 
at the impurity level. The critical temperature in the 
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FIG. 2: Phase diagrams for the 'activated' RKKY-BMP per- 
colation combined model in two typical cases. The dashed 
lines represent the two critical temperatures of the 'activated' 
RKKY model, Eq. ((2)|. The solid line is the critical tem- 
perature in the bound magnetic polaron percolation model 
given by Eq. The shaded areas are the ferromagnetic 

phases and the white ones are paramagnetic. The parame- 
ters used are rii = 3 x 10^^ cm~"^, J = 1 eV, ao — 3.23A 

.19 „™-3 ,X (a|nco ~ 0.002), (b) 

0.0125). 



and (a) Uco = 3 x 10"" cm ", as 
Uco = 1 X W^^ cm~^, as = 5A (a| 



rico 



BMP percolation model is given by [ic 



CLB 



(3) 



where as is the carrier localization radius, and n* is the 
carrier density at Tf""''^. Eq. ([3]) is valid in the low car- 
rier density limit a^n* << 1. As the density of carriers 
involved in polaron formation n* increases with A, T'^'^^^ 
also increases, saturating for large values of A, when 
~ "^cO; ^ shown in Fig. [5] (solid line). The value of 
jiporc jj^a^jj^iy dominated by the value of a^n* , on which 
it depends exponentially. Another mechanism that could 
arise when there are no free carriers is the Bloembergen- 
Rowland-type mechanism proposed in Ref. [l6| in the 
context of Mn-doped III-V semiconductors, where the 
virtual electron excitations take place between the im- 
purity levels and the band. We do not expect, however, 
this mechanism to be stronger than the BMP percolation, 
and therefore, it does not affect our conclusions. 

In Fig. [2] we show two typical phase diagrams resulting 
from the interplay of the 'activated' RKKY and the BMP 
percolation models. The parameters chosen here are rep- 
resentative of real systems. We have fixed the magnetic 
impurity concentration to a; = 0.1, J = 1 eV, and 5* = 1. 




FIG. 3: The shaded area gives the values of as and ricO 
compatible with the re-entrant EM behavior illustrated in 
Fig. (Ha). The curve is given by a%nco = 0.00824. 



The other two main parameters left are n^o and as ■ rico is 
chosen so that nco/n-i < 0.1, as it is usually found in DMS 
systems due to strong carrier compensation. Finally, as 
is chosen to fulfill the applicability conditions of Eq. 
There are two qualitatively different phase diagrams that 
can occur depending on the value of a^yico- For the low- 
est values of o^tIcq, T3 = rP°''^ is relatively small, and 
we get the novel scenario depicted in Fig. [2] (a): in region 
I, corresponding to A/T^^^^ ^ 0.32 for the parameter 
values chosen in this figure, the system is FM with critical 
temperature Ti; in region III, with A/T^^^ > 0.73 for 
any set of parameters, the system is ferromagnetic with 
Tc ^ T3] finally, in region II, which corresponds to inter- 
mediate values of A (0.32 < A/T5,kky < 0.73), T3 < T2 
and the system shows re-entrant FM with three distinct 
critical temperatures Ti, T2, and T3. For larger values 
of a%nco (but stiU << 1), as in Fig. [H (b), T3 > T2 
for all values of A and the system is ferromagnetic for 
T<max(ri,T3). 

Although our theoretical analysis, being physically mo- 
tivated, is formally correct, the question naturally arises 
about the observability of our predicted re-entrant FM 
behavior in DMS materials. There are three general con- 
ditions that should be met: (i) at T = there should 
be no carriers in the conduction (or valence) band so all 
the carriers mediating RKKY come from thermal activa- 
tion; (ii) for the scenario illustrated in Fig. [2] (a) to occur 
we need to fulfill the condition TP"'=(nco) < T^'""" ~ 
0.3 X T^^^^(nco); and (iii) we also have to ensure that 
the system is in region II in Fig. [5] (a) to actually observe 
re-entrant FM. 

Condition (i) basically implies that the system has to 
have activated-like resistivity (decreasing as T increases) 
limiting the suitable systems to those with so-called 
'insulating behavior', for example, very lightly doped 
(Ga,Mn)As, (In,Mn)As, and magnetic oxides. Condition 
(ii) puts restrictions on the relative values of and rico- 
The equation in (ii) is independent of J and rii and re- 
duces to the inequality a^Uco < 0.00824 for 5 = 1. This 
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result is illustrated in Fig. [3] where the curve is given by 
the maximum allowed value of the product a^rico. Be- 
low this curve we get the scenario in Fig. [2] (a) while the 
scenario in Fig.[2](b) occurs for gb and rico values above 
the curve. The carrier confinement radius is given by 
ub — e{m/m*)a where a — 0.52 A, e is the semiconductor 
dielectric constant, and m* is the polaron mass (usually 
m*/m > 1). For diluted magnetic semiconductors and 
magnetic oxides we generally expect 3A< ub < 10 A. 
Hence, as shown in Fig. [3l which is the important ma- 
terials phase diagram to keep in mind in searching for a 
suitable system to observe our predicted re-entrant FM, 
we need to keep a relatively low density of carriers, but 
well within experimentally achievable values. 

Condition (iii) can be fulfilled, in general terms, when 
A is comparable to Tc- This rules out Gai_a;Mna;As, with 
A/Tc ~ 10 (Tc - 100 K, and A ~ 110 meV [I^), which 
would rather be in region III (for very low carrier density 
samples). On the other hand, the re-entrant FM scenario 
might be applicable in Ini-ajMn^jAs, which has very low 
activation energies (see, for example. Fig. 2 in Ref. 
and low Tc ^ 50 K. However, the most likely candidate 
for observing our predicted re-entrant FM is possibly a 
doped magnetic oxide material, e.g. Tii_a:Coa;02 (with 
30meV < A < 70 meV [19] and very similar critical tem- 
peratures Tc ~ 700 K), which is an insulator at low T, 
but exhibits essentially T-independent Drude transport 
behavior, due to thermally activated carriers, at room 
temperature or above. In such a system, it is quite pos- 
sible that the high-T FM is mediated by thermally acti- 
vated carriers, whereas at low T a BMP percolation FM 
takes over [20| as in region I of Fig. [2] (a). 

We have neglected the direct exchange interaction be- 
tween the magnetic impurities, which is short-ranged, 
and could also play a role, particularly when the den- 
sity of impurities rij is large, completely destroying the 
lower FM phase 2l|, and leading to a spin-glass phase. 
These spin-glass low-T phases have been observed in II- 
VI DMS [3, and (Ga,Mn)N [H, in samples that never 
show FM order. More relevant to our case, a suppre- 
sion of magnetization at low T in zero field cooled curves 
has been reported in the magnetic oxide V-doped ZnO 
and possibly in other magnetic oxides ^23|. In princi- 
ple, therefore, direct antiferromagnetic exchange between 
the magnetic dopant impurities could compete with (or 
even suppress) our predicted low-T re-entrant FM phase, 
but the origin of this direct exchange being completely 
different from the carrier-mediated mechanisms produc- 
ing the re-entrant FM behavior itself, we find it difficult 
to believe that such a suppression of re-entrance can be 
generic. We suggest detailed T-dependent magnetization 
studies in the shaded region of our Fig. [3] to search for 
our predicted re-entrant FM behavior. 

We have considered, using physically motivated ef- 
fective models of FM in DMS materials, the intrigu- 
ing possibility of generic re-entrant FM in an insulat- 



ing class of DMS systems. The models we use are con- 
sidered to be successful minimal carrier-mediated mod- 
els 0, S, S, [13, [m for understanding FM behavior and 
predicting ferromagnetic Tq in itinerant and localized 
DMS materials. The new idea in this work has been 
to point out that these 'competing' FM mechanisms, e.g. 
'activated' RKKY and BMP percolation, could, in prin- 
ciple, exist together in a single sample where thermal 
activation leads to a high-T effective free-carrier mecha- 
nism mediating the RKKY ferromagnetism, and the low- 
T FM, where thermal activation of free carriers is expo- 
nentially suppressed, is mediated by localized bound car- 
riers in an impurity band through the polaron percolation 
mechanism. We show that, depending on the materials 
parameters, such a situation with competing high-T and 
low-T FM mechanisms generically allows for re-entrant 
FM with an intermediate-temperature anomalous para- 
magnetic phase intervening between the higher tempera- 
ture RKKY FM phase and the lower temperature BMP 
percolation FM phase. More experimental work in DMS 
materials is needed to confirm the existence of re-entrant 
FM, but the phenomenon should exist on firm theoretical 
grounds. 

This work is supported by the NSF and the NRI SWAN 
program. 
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